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The solution to this problem can be obtained by solving the linearized unsteady Euler equations. Let 
the unsteady flow field be given by 


U (x, t) = Up(x) + u(x, t ) 

(1) 

p(x,t) = p 0 (x) +p'(x,t) 

(2) 

p(x,t) = p 0 (x) + p\x,t) 

(3) 

s(x, t) = s 0 + s'(x, t) 

(4) 


where the entropy s 0 is constant, and u, p ' , //, and s' are the unsteady perturbation velocity, pressure, 
density and entropy, respectively. Zero subscripts denote mean flow quantities which are assumed to be 
known. 

Substituting (1) - (4) into the nonlinear Euler equations and neglecting products of small quantities, 
one obtains the linearized equations 
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where ^ + U 0 • V is the convective derivative associated with the mean flow. 

If the mean velocity U 0 can be expressed as the gradient of a potential then equations (5) - (7) can 
be reduced to a single, non-constant coefficient, inhomogeneous convective wave equation (refs. 1,2) 
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where the unsteady velocity is decomposed into a known vortical component and an unknown potential 
component V0, 

u(x,t) = - bV</>. (9) 

The unsteady pressure is given by 

p = ( 10 ) 

An unsteady aerodynamic code, called GUST3D (ref. 3), has been developed to solve equation (8) for 
flows with periodic vortical disturbances. The code uses a frequency-domain approach with second-order 
central differences and a pressure radiation condition in the far field. GUST3D requires as input certain 
mean flow quantities which are calculated separately by a potential flow solver. This solver calculates the 
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mean flow using a Gothert’s Rule approximation (ref. 3). On the airfoil surface, it uses the solution calculated 
by the potential code FL036 (ref. 4). Figures 1-2 show the mean pressure along the airfoil surface for the 
two airfoil geometries. 

In Figures 3 - 8, we present the RMS pressure on the airfoil surface. Each figure shows three GUST3D 
solutions (calculated on grids with different far- field boundary locations). Three solutions are shown to 
provide some indication of the numerical uncertainty in the results. 

Figures 9-13 present the acoustic intensity. We again show three solutions per case. Note that no 
results are presented for the Aq = = 2.0 loaded airfoil case, as an acceptable solution could not be obtained. 

A few comments need to be made about the results shown. 

First, since the last Workshop, the GUST3D code has been substantially upgraded. This includes imple- 
menting a more accurate far-held boundary condition (ref. 5) and developing improved gridding capabilities. 
This is the reason for any differences that may exist between the present results and results from the last 
Workshop. 

Second, the intensity results on the circle R = 4(7 were obtained using a Kirchoff method (ref. 6). The 
Kirchoff surface was the circle R = 2 C. 

Finally, the GUST3D code is most accurate for low reduced frequencies. A new domain decomposition 
approach (ref. 7) has been developed to improve accuracy. Both the single domain and domain decomposition 
approaches were used in generating the present results. 
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|p| A 2 sin(0) |p| A 2 sin(0) 



Figure 9. a Acoustic intensity on circle R = 1 C, 
Case 1, kl = k2 = 0.1. 



Figure 9.b Acoustic intensity on circle R = 2 C, 
Case 1, k 1 = k2 = 0.1. 
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|p| A 2 sin(0) |p| A 2 sin(0) 



Case 1, kl = k2 = 0.1. 



Figure 10. a Acoustic intensity on circle R = 1 C, 
Case 1, kl = k2 = 1.0. 
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|p| A 2 sin(0) |p| A 2 sin(0) 



Figure lO.b Acoustic intensity on circle R = 2 C, 
Case l,kl =k2 = 1.0. 



Figure lO.c Acoustic intensity on circle R = 4C, 
Case l,kl =k2= 1.0. 
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|p| A 2sin(0) |p| A 2sin(0) 



Figure 1 l.a Acoustic intensity on circle R = 1 C, 
Case 1, kl = k2 = 2.0. 



Figure 1 1 .b Acoustic intensity on circle R = 2 C, 
Case 1, kl = k2 = 2.0. 
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|p| A 2 sin(0) |p| A 2 sin(0) 



Figure 1 l.c Acoustic intensity on circle R = 4C, 
Case 1, kl = k2 = 2.0. 



Figure 12.a Acoustic intensity on circle R = 1 C, 
Case 2, kl = k2 = 0.1. 
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|p| A 2 sin(0) |p| A 2 sin(0) 



Figure 12.b Acoustic intensity on circle R = 2 C, 
Case 2, kl =k2 = 0.1. 



Figure 12.c Acoustic intensity on circle R = 4 C, 
Case 2, kl = k2 = 0.1. 
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|p| A 2 sin(6) |p| A 2 sin(0) 



|p| A 2 sin(O) 

Figure 13. a Acoustic intensity on circle R = 1 C, 
Case 2, kl = k2 = 1.0. 



Figure 13.b Acoustic intensity on circle R = 2 C, 
Case 2, k 1 = k2 = 1.0. 
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|p| A 2 sin(6) 



Figure 13.c Acoustic intensity on circle R = 4 C, 
Case 2, kl = k2 = 1.0. 
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